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ABSTRACT 

We present broad-band and spectropolarimetry of the Type la SN 1996X obtained 
on April 14, 1996 (UT), and broad-band polarimetry of SN 1996X on May 22, 1996 
(UT), when the supernova was about a week before and 4 weeks after optical maximum, 
respectively. The Stokes parameters derived from the broad-band polarimetry are con- 
sistent with zero polarization. The spectropolarimetry, however, shows broad spectral 
features which are due intrinsically to an asymmetric supernova atmosphere. The spec- 
tral features in the flux spectrum and the polarization spectrum show correlations in 
the wavelength range from 4900 A up to 5500 A. The degree of this intrinsic component 
is low (~ 0.3%). Theoretical polarization spectra have been calculated. It is shown 
that the polarization spectra are governed by line blending. Consequently, for similar 
geometrical distortions, the residual polarization is smaller by about a factor of 2 to 
3 compared to the less blended Type II atmosphere, making it intrinsically harder to 
detect asphericities in SNIa. Comparison with theoretical model polarization spectra 
shows a resemblance to the observations. Taken literally, this implies an asphericity of 
~ 11 % in the chemical distribution in the region of partial burning. This may not 
imperil the use of Type la supernovae as standard candles for distance determination, 
but nontheless poses a source of uncertainty. SN 1996X is the first Type la supernova 
for which spectropolarimetry revealed a polarized component intrinsic to the supernova 
and the first Type la with spectropolarimetry well prior to optical maximum. 



Subject headings: stars: individual (SN 1996X) - stars: supernovae - stars: polarimetry 



^lso Beijing Astronomical Observatory, Beijing 100080, P. R. China 

2 email : lifan® astro . as . utexas . edu 

3 email : wheel® astro . as . utexas . edu 
4 email: pah@alla.as.utexas.edu 

5 Also: Institute for Theoretical Physics, University of Basel, Klingelbergstr. 82, CH-4056 Basel, Switzerland 



1 



1. Introduction 

Polarimetry is a powerful way of probing asym- 
metries in supernova explosions. Because polariza- 
tion produced within the supernova atmosphere may 
have spectral features that are distinguishable from 
the featureless interstellar polarization (ISP), which 
is caused by dust within both the Galaxy and the host 
galaxy spectropolarimetry is usually a more power- 
ful tool than broad-band polarimetry in decomposing 
the intrinsic polarization of a distant supernova from 
ISP. 

SN 1987A and SN 1993J are the only two super- 
novae for which spectropolarimetry has been pub- 
lished and intrinsic polarizations are positively de- 
tected. SN 1987A is so far the best observed su- 
pernova for which linear polarimetry was obtained 
from several days to about 260 days after explo- 
sion (Mendez et al. 1988; Cropper et al. 1988; 
Jcffcry 1991a). The degree of polarization evolved 
with time, indicating that the cause of the polar- 
ization is related intrinsically to SN 1987A. The 
data were analyzed (Hoflich 1987; Jeffery 1991b) in 
terms of the photospheric scattering model (Brown 
& McLean 1977; Shapiro & Sutherland 1982). Re- 
cently, Wang & Wheeler (1996a) provided a differ- 
ent view in which time delayed scattering by a hy- 
pothesized circumstellar dust clump successfully re- 
produced both the broad-band and spectropolarime- 
try of the supernova and the early infrared light curve 
of SN 1987A (Bouchet et al. 1989). Linear polariza- 
tion indicated by polarization changes across spectral 
features was also detected in SN 1993J (Trammell, 
Hines, & Wheeler 1993). Broad band polarimetry 
shows also variable polarization before and after the 
second optical maximum of SN 1993J (Doroshenko, 
Efimov, & Shakhovskoi 1995). 

Only two Type la supcrnovae - SN 1983G (McCall 
et al. 1984) and SN 1992A (Spyromilio & Bailey 1992) 
have been previously observed by spectropolarimetry. 
The noise levels of the SN 1983G and SN 1992A data 
are 0.5% and 0.3%, respectively. These data indicate 
SN 1983G and SN 1992A arc not polarized at a level 
higher than 0.5%. 

As part of our program of systematic supernova po- 
larimetry, we have observed 5 supcrnovae with broad- 
band polarimetry and compiled a catalog of all the su- 
pcrnovae with polarimetry reported in the literature 
(Wang et al. 1996). This sample shows that all the 
Type II supernovae with sufficient data are intrinsi- 



cally polarized while no intrinsic polarization can be 
established for any of the Type la supernovae. 

Supernovae, of both Type la and II, are now being 
widely used as standard candles for distance determi- 
nations. As a self-consistency test, it is imperative 
to set some observational constraints on the degree of 
the polarization. Furthermore, normal spectroscopy 
provides only part of the information carried by the 
supernova light. Spectropolarimetry can further con- 
strain various models for spectral line formation and 
radiative transfer through supernova atmospheres, as 
shown in recent examples by Hoflich (1995a), Hoflich 
et al. (1996), and Hoflich, Wheeler, & Wang (1997). 
Furthermore, it may be the only way to detect chem- 
ical blobs in the ejecta which usually produce little 
observable effects on the flux spectra of a supernova 
(Hoflich, Wheeler, & Wang 1997). 

2. The Observations 

SN 1996X in NGC 5061 was discovered indepen- 
dently by several observers (Garradd 1996) on April 
12.5 UT at a V magnitude about 13.5. Spectroscopy 
on April 14 UT showed SN 1996X to be a Type la 
supernova prior to optical maximum (Suntzeff 1996; 
Wang and Wheeler 1996b; Benetti and Patat 1996). 
Our spectropolarimetry of SN 1996X was obtained 
on April 14.3 using the Imaging Grism Polarimctcr 
(IGP) mounted at the Cassegrain focus of the 2.1 me- 
ter telescope of the McDonald Observatory. The IGP 
is a simple, high efficiency dual beam polarimeter 
which can be easily switched between imaging and 
spectropolarimetry mode. A rotatable half waveplate 
was used as the polarization analyzer. A description 
of the instrument can be found in Hill & Trammell 
(1996) and Trammell (1994). A TK4 (1024 x 1024) 
CCD was used as the detector. A broad filter with 
pass band from 4400 to 7500 was used for the broad- 
band polarimetry. Such a filter effectively filters out 
photons with wavelengths beyond the effective range 
of the waveplate, but allows a large number of pho- 
tons to be transmitted for polarimetry. The slit width 
used for the spectropolarimetry was 2".l which gave 
a spectral resolution of about 14A. Wavelength cali- 
bration was obtained by taking exposures of an argon 
lamp. 

The data were taken with clear and dark sky. Sev- 
eral unpolarized and polarized standard stars were 
observed each night. The data reduction process is 
the same as outlined in Miller, Robinson, & Goodrich 
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Fig. 1. — Instrumental polarization from observa- 
tions of unpolarized standard stars. 
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Fig. 2.— (a) The flux spectra of SN 1996X on April 
14.3, 1996 UT in arbitrary units; (b) theNStokes pa- 
rameter Q; (c) the Stokes parameter U; (d) the noise 
level of the Stokes parameters. In (b) - (c), the thin 
lines are the observed data at the original sampling 
step of 3.8A/pixel, the thick lines are the same data 
resampled with a sample window width of 64.8A(17 
pixels). 




Fig. 3. — The wavelength range showing correlations 
between the flux spectrum (dotted line), the polar- 
ization spectrum from the resampled data (thin solid 
line), and the polarization of the data smoothed by 
convolving with a Gaussian of FWHM of 64.8A(thick 
solid line). 



(1988). The instrumental polarization obtained by 
observing several unpolarized standards is shown in 
Figure 1 (a) and 1 (b) and is found to be less than 
0.1% in the entire wavelength range. The polarized 
standards were used to calibrate the polarization an- 
gles. In addition, observations were also obtained by 
inserting a polarizer at the top of the waveplate to 
characterize the efficiency of the waveplate and the 
wavelength dependence of the fast axis of the wave- 
plate. The polarization efficiency is typically 95% at 
around 6500 A and decreases to about 90% at 4500 
A. The variation of the fast axis with wavelength is 
less than 2°. These effects and also the instrumental 
polarization were taken into account and corrected for 
in the data that will be presented below. 

3. Data and Results 

Although the exact phase of the supernova awaits 
the publication of an accurate light curve, the super- 
nova was at pre-maximum phase on April 14, as indi- 
cated by the flux spectrum shown in Figure 2 (a). The 
continuum was very blue and the Fe II lines were still 
not fully developed. The Stokes parameters Q and 
U derived from the spectropolarimetry are shown in 
Figures 2 (b) and 2 (c), and their errors in Figure 2 
(d) . The dominant error in Q and U was from photon 
statistics, and was the largest in the red portion of the 
spectra where it reached 1%; typical errors were less 
than 0.6% in regions around 5500 A. 

Polarizations higher than 1% can be reliably ex- 
cluded from the Q and U spectra shown in Figure 2. 
However, at a lower level of about 0.3%, some broad 
features are present. For a simple test, the Q and U 
spectra were convolved with a Gaussian of FWHM 
of 65 A. As shown in Figure 3, in the wavelength 
range from 4900 to 5500 A, there exists a correlation 
between the wavelengths of these broad features and 
the spectral features in the flux spectrum. 

To investigate the authenticity of this polarization 
variation with wavelength, the data were resampled 
to increase the signal to noise ratio. This was done 
by deriving the Stokes parameters within a certain 
wavelength window AA centered on a particular wave- 
length, A. This is equivalent to calculating the follow- 
ing weighted mean of the original Stokes parameters: 

r-X+AX/2 a+AA/2 

Q(A|AA) = / N{X')Q{X-X')d\' / / N(X')dX', 

JX-AX/2 JX-AX/2 

(1) 

where N is the number of counts per wavelength 
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interval. The corresponding equation for U can be 
similarly derived. Such an operation increases the 
signal to noise ratio of the data, but at the expenses 
of degrading the spectral resolution. Fortunately, the 
spectral features in a supernova spectrum are in gen- 
eral as broad as 150 A, and useful information can 
still be obtained even if the resolution is as low as 100 
A. The Stokes parameters with an integration win- 
dow of width 2 x AA = 64.8 A are shown in Figure 2 
(b) and (c). The noise level (a) which is also shown 
in Figure 2 (d), is now typically around 0.1%. 

The degree of polarization P can be constructed 
from the Stokes parameters Q and U . However, it is 
well known that for relatively high noise level, P = 
\lQ 2 + U 2 is biased toward values larger than the 
true degree of polarization P Q ; P is usually not the 
best estimate of P . The distribution function for 
(P, 9) takes the form (Simmons & Stewart 1985) 

lira la a a a 

(2) 

where (P ,9 ) stands for the true degree of polar- 
ization and polarization position angle. Simmons & 
Stewart (1985) constructed optimal estimates of de- 
gree and position angle of polarization (P and 9 D ) 
from the marginal distribution of polarization ob- 
tained by integrating over 9 in equation (2). This 
method, however, does not usually simultaneously 
optimize P Q and 9 Q . A simpler method can be con- 
structed by requiring the best estimate of (P ,9 ) to 
take the values for which (P, 9) gives the maximum 
of the distribution f(P,9). It is easy to derive from 
equation (2) that this results in the following simple 
formula for P and 9 a 

P = P-(J 2 /P P > a 

P = P < a (3) 

8 = 9 

Both the polarization P a derived from equation (3) 
and P are shown in Figure 4(a). They differ mainly 
in wavelength regions with low signal to noise ratios. 
The polarization position angles could not be reliably 
determined for data with large errors, therefore only 
position angles calculated from the resamplcd Q and 
U are shown (Figure 4(b)). 

An enlarged version of the part of the spectrum 
with strong correlations between the spectral features 
in the polarimetry and flux spectrum is shown in Fig- 
ure 3, where the features are marked by dashed lines 
for illustrative purposes. The corresponding spectral 



features in the polarization spectrum which are cor- 
related with the spectral lines in the flux spectrum 
are blue shifted by about 2200 kms -1 . No strong 
correlation between the flux and polarization spec- 
tra could be established at wavelengths longer than 
6000 A. This may be due to some intrinsic character- 
istics of the polarization of Type la SN (see §4 for a 
more detailed discussion). The position angles change 
across spectral features. This should not be a surprise 
when compared with the best observed supernova, SN 
1987A, where the intrinsic component shows dramatic 
variations of both degree and polarization position an- 
gle across spectral features (Cropper et al. 1988). 




5000 6000 7000 

Wavelength (A) 



Fig. 4. — (a) The number of counts in a single beam 
of a 20 minute exposure on a logarithmic scale (dot- 
ted line); the degree of polarization resampled with a 
sample window width of 64. 8A; the degree of polar- 
ization P Q as defined in equation (2) (thick solid line), 
(b) The polarization position angle of the resampled 
data. 

It is impossible to accurately correct for the effect 
of ISP. The Galactic extinction to the host galaxy 
NGC 5061 is Ab = 0.25, which alone can pro- 
duce polarizations as large as 0.7% according to the 
limit derived by Serkowski (1970). The ISP correc- 
tion is even more complicated by the polarization due 
to dust in the host galaxy; however, the ISP follows 
Serkowski's law, and is a smooth function of wave- 
length; the uncertainties introduced by the ISP will 
not change the Q and U spectral structures on wave- 
length scales less than a few hundred A. The ISP 
correction can, however, dramatically change the ap- 
pearance of the polarization spectra. The spectral 
features in Figure 4(a) may change from 'emission' to 
'absorption' even if a very small amount of ISP cor- 
rection ~ 0.2% is required. This should be borne in 
mind when interpretating the polarization spectra. 

In addition to the spectropolarimetry discussed 
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above, we have also obtained broad-band polarimetry 
using a wide filter with pass band from 4000 - 7500 
A. Three and five data sets were obtained on April 
14 and May 22, respectively. Each data set consists 
of four 100 second exposures which gave estimates of 
the Stokes parameters. The measured Q and U are 
practically zero for all data sets. The averaged values 
are Q = 0.079 ± 0.041% and U = 0.056 ± 0.039% 
on April 14, and Q = 0.090 ± 0.021% and 
U = 0.063 ± 0.019% on May 22, where the errors are 
due to photon statistics. The systematic instrumen- 
tal uncertainty is around 0.06% (as indicated by re- 
peated observations of unpolarizcd standards) and is 
therefore the dominant error for the broad-band data. 
The broad-band data are consistent with zero polar- 
ization. We do not confirm the short-term temporal 
variations of polarization reported in early polarime- 
try of another Type la SN 1972E (Lee, Wamsteker, 
& Wisniewski 1972, quoting observations made by 
Serkowski and Wamsteker). 

It should be pointed out that the null detection 
in broad-band polarimetry did not conflict with the 
spectral features detected in the spectropolarimetry 
data. As a check for consistency, we have weighted the 
spectropolarimetry data with the transmission curve 
of the filter and the photon counts at different wave- 
lengths to simulate broad-band polarimetry. The re- 
sulting Q and U arc in agreement with the broad-band 
polarimetry data. 

4. Models 

Besides the attempt to identify features of the po- 
larization with those of the flux spectrum, a more di- 
rect approach is a comparison with theory which pro- 
vides additional information on the spectral patterns 
to be expected. A simplified approach is suitable to 
answer the following questions: what does the polar- 
ization spectrum look like for a SN la, what aspheric- 
ity is needed to be consistent with the observation of 
SN1996X, and how do the predicted spectra compare 
with the observations ? Delayed detonation models 
(Khokhlov 1991) have been found to provide a good 
representation of the observed light curves and spec- 
tra of normal bright SN la (Hoflich 1995b, Hoflich 
& Khokhlov 1996 and references therein) and we use 
such a model to produce a representative theoreti- 
cal SN la polarization. The polarization spectrum 
is given by Monte Carlo calculations based on the 
following assumptions: a) homologously expanding, 



ellipsoidal envelopes with density and chemical pro- 
files given by a hydrodynamical model, b) occupation 
numbers given by local thermodynamical equilibrium 
(LTE); c) electron scattering, bound- free and free- free 
for continuum opacities; d) lines treated in a Sobolev 
approximation with an assumed constant thermaliza- 
tion fraction; e) line transitions result in depolariza- 
tion; f) the temperature structure is given; g) and 
elliptical geometry with a axis ratio being indepen- 
dent of radius. Note that discrepancies between pre- 
dicted and observed wavelengths of polarization fea- 
tures corresponding to 1000 to 2000 kms -1 may be 
expected because this difference is well within the ob- 
served variation of expansion velocities among normal 
bright SNe la. Lines of the same ion must, neverthe- 
less, be expected to be shifted consistently. For more 
details, see Hoflich et al. (1995). 



SN1996X vs. Model DD200c at day 10 
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Fig. 5. — Comparison of the theoretical polarization 
spectrum of the delayed detonatinon model DD200c 
five days before maximum light (thick line; Hoflich et 
al. 1996; axis ratio 0.88 seen at an inclination of 30°) 
with the observations (thin line). The identification of 
some features are given. Most of the weaker features 
are due to iron group elements in the second and third 
ionization stage. The S II and Si II features at about 
4900, 5300, 5600 and 5670 A are too strong and the 
overall polarization above 6500 A is too large. Possi- 
ble solutions are discussed in the text. Note, that the 
observed spectrum is oversampled although smoothed 
with a triangle of 64.8 AFWHM. 

For the density and chemical structure we used 
model DD200c with T(r) based on light curve calcula- 
tions corresponding to about 5 days before maximum 
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(Hoflich et al. 1996). In Figure 5, the polarization 
spectrum for an oblate ellipsoid with an axis ratio of 
0.89 is given which is seen at an inclination of 30° 
which fits best. Note that, quantitatively, the spec- 
trum depends critically on the inclination. The depo- 
larization in the ubiquitous metal lines has two main 
effects. Firstly, the mean polarization is about a fac- 
tor of 2 to 3 smaller compared to that of a Thomson 
scattering atmosphere which provides a good approx- 
imation to Type II supcrnovae. Consequently, the 
same asphericity will produce much less polarization 
for an SN la compared to an SN II. Secondly, P varies 
strongly with wavelength on a typical scale of about 
the Doppler width. The variations with wavelength 
are significantly stronger than those in the flux spec- 
tra because it takes only one interaction to depolarize 
but many to thcrmalizc a photon. 

A comparison with the observations shows some 
qualitative similarity of the wavelength distribution 
of the features. The pattern in the models is mainly 
influenced by the atomic physics rather than the de- 
tailed explosion model. Therefore, the agreement or 
disagreement may be used to ask if the measurement 
can be regarded as a real detection. Most of the ob- 
served features have their equivalents in the theoreti- 
cal model, although the predicted polarization above 
6500 A is much higher than observed (see below). The 
mean size of P implies an asphericity of « 1 1 % at an 
inclination angle of 30°for the symmetry axis of the 
ejecta (Figure 5). 

There are also some distinct shortcomings in the 
model polarization spectrum. The features corre- 
sponding to the products of partial burning (S II and 
Si II) are somewhat too strong compared to the ob- 
servations. In principle, this may be compensated by 
slightly higher temperatures, small changes in the ex- 
plosion model (e.g. a slightly larger Ni production) 
or, as test calculations have shown, by making the 
outer Si-rich layers almost spherical. A more severe 
problem is the spectrum above 6500 A. Although the 
spectral features resemble the models, the absolute 
polarization is much lower. A component due to inter- 
stellar polarizaton (U=-0.1 %, Q = 0.25 %) helps to 
bring up P above 6500 % but completely destroys any 
reasonable fit at other wavelengths and, thus, does 
not solve the problem. This ISP is also inconsistent 
with broad-band data (see above). The increasing 
predicted polarization is insensitive to the details of 
the explosion models. It is caused primarily by the 
decreasing line blending towards longer wavelengths. 



Lower line blending also implies that the polariza- 
tion in the red is produced at deeper layers (Hoflich 
1995b). Consequently, a possible solution to the prob- 
lem is a nearly spherical, inner Ni region correspond- 
ing to expansion velocities w 8000km s -1 . This, to- 
gether with the depolarization in Si and S, imply that 
asphericities must be mainly attributed to the distri- 
bution of chemical elements in the transition region 
between the Ni and Si layers (Hoflich et al. 1996). 

The correlation between flux and polarization spec- 
tra depends sensitively on the size of the line blan- 
keting, thermalization parameter of lines and the ac- 
tual wavelength. For a more detailed discussion see 
Hoflich, Wheeler, & Wang (1997). 

5. Discussions and Conclusions 

SN 1996X is the first Type la supernova for which 
spectropolarimetry has suggested a polarized compo- 
nent intrinsic to the supernova. The peak to valley 
spectral variation in the Q and U spectra (Figure 4) 
is as high as 0.6% in many cases which, at an error 
level of 0.15%, gives a 4cr detection of the spectral 
features. The degree of the intrinsic polarization can 
be measured in terms of the deviation of the spectral 
features from the wavelength averaged degree of po- 
larization. For SN 1996X, this is about 0.3%, which 
is small compared with the typical degree of polar- 
ization in SN 1987A and other Type II SNe (Wang 
et al. 1996). Previous observations of SN 1983N and 
SN 1992 A have given only upper limits on the degree 
of polarization of Type la supernova (Spyromilio & 
Bailey, 1992; McCall et al. 1984). The SN 1996X 
observations were obtained prior to the optical max- 
imum, and the detected polarization is most likely to 
be due to a distorted photosphere or element distribu- 
tion. The correlation of features between the polar- 
ization spectra and flux spectra is not perfect in the 
entire observed wavelength range. This is expected 
because such a correlation depends sensitively on mul- 
tiple scattering effects, the ratio between continuum 
and line optical depths, and thermalization in lines. 
The 2200 kms -1 blue shift of the features in the po- 
larization spectrum relative to the flux spectrum (cf. 
Figure 3) can also be understood as due to multiple 
scattering. The model calculation shows that a strong 
feature in the polarization spectra does not require 
the presence of a strong feature in the flux spectra. 
Lines which hardly change the flux may produce sig- 
nificant depolarizations. The comparison with theo- 
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retical polarization spectra suggests an asphericity of 
about 11 % in the distribution of chemical elements. 
There are evidences that, the strong asphericity is 
limited to the transition region between nuclear sta- 
tistical equillibrium (NSE) and partial burning. The 
distorted nature of the supernova envelope might arise 
from inhomogenous burning during the deflagration 
phase of the burning front (Khokhlov 1995). 

Such a low degree of polarization probably does not 
seriously endanger the use of Type la as calibrated 
candles for distance indicators, but nonetheless poses 
a source of uncertainty. We will discuss these effects 
including the effects of inclination and the relation 
between flux and polarization spectra quantitatively 
in a separate study (Hoflich, Wheeler, & Wang 1997). 
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